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Abstract. In Schaerer et al. (1995, Paper I) we have pre- 
^e^ited the first "combined stellar structure and atmo- 
^^here models" (CoStar) for massive stars, which consis- 
tently treat the entire mass loosing star from the center 
T oTrt to the outer region of the stellar wind. The models use 
(UQto-date input physics and state-of-the-art techniques to 
Sbdel both the stellar interior and the spherically expand- 
•^jE^ non-LTE atmosphere. The atmosphere models include 
^jne blanketing for all elements from hydrogen to zinc. 

The present publication covers the spectral evolution 
corresponding to the main sequence interior evolution dis- 
CSussed in Paper I. The CoStar results presented in this 
CjQper comprise: (a) flux distributions, from the EUV to 
le far IR, and the ionizing fluxes in the hydrogen and 
hum continua, (b) absolute optical and infrared UB- 
IIJHKLMN photometric magnitudes and UV colors, (c) 
©^tailed line blanketed UV spectra, and (d) non-LTE hy- 
"^jb gen and helium line spectra in the optical and IR, in- 
ghid ing theoretical K band spectra. These results may, 
eig., be used for population synthesis models intended to 
lyidy the massive star content in young starforming re- 
"l^pns. 

^ We compare our results with other predictions from 
LTE and non-LTE plane parallel models and point out the 
improvements and the importance of using adequate at- 
mosphere models including stellar winds for massive stars. 
Particular emphasis is given to comparisons of the UV 
spectral evolution with observations, including continuum 
indices and several metal line signatures of P-Cygni lines 
and broad absorption features. Good agreement is found 
for most UV features. In particular, we are able to repro- 
duce the strong observed Fe III 1920 A feature in late O 
and early B giants and supergiants. This feature is found 



Send offprint requests to: D. Schaerer 



to depend sensitively on temperature and may be used to 
derive effective temperatures for these stars. 

We also derive a simple formula to determine mass loss 
rates from the equivalent width of hydrogen recombination 
lines (Ha, Pa and Ba) for OB stars showing net emission 
in one or more of these lines. 

Key words: Stars: atmospheres - early-type - evolution 
- fundamental parameters - Hertzsprung-Russel (HR) di- 
agram - mass-loss 



1. Introduction 

This is the second paper on "combined stellar struc- 
ture and atmosphere models" (CoStar). In Co Star models 
we treat the stellar interior, the line-blanketed non-LTE 
spherically extended atmosphere and the wind in its en- 
tirety. The models are intended to study massive stars 
(M 2; 2OM ) as for these stars mass loss through a stel- 
lar wind is important for the stellar evolution (cf. Chiosi 
& Maeder 1986, Maeder & Conti 1994) and also because 
of the impact of the winds of these stars on the emergent 
flux distribution (e.g. Kudritzki & Hummer 1991). CoStar 
models may therefore be used to study the continuum and 
(metal)line spectra of stars in evolutionary phases corre- 
sponding to OB, LBV, Of/WN, and Wolf-Rayet (WR) 
stars. 

An important advantage and motive for this work is 
that our approach allows for a more 'direct' comparison 
of stellar evolution results with observations. A second 
incentive for this work is that our method allows us to 
study the effect of a spherically expanding atmosphere 
and stellar wind on the envelope, interior structure and 
evolution of massive stars. Last, but not least, our models 



e.g. WR winds and LBV outbursts, which are likely to 
require a consistent modeling of both the stellar envelope 
and the wind. 

Our approach implies the consistent modeling of both 
the UV, optical and IR spectral ranges. A future aim is to 
derive stellar parameters from those diagnostics that are 
most sensitive to a specific parameter. This includes the 
determination of metal abundances from UV spectral line 
features. These UV features are at the same time likely to 
yield valuable constraints on the metal ionization, which 
is important for the dynamical modeling of stellar winds. 

In the present work we concentrate on the main se- 
quence (MS) and present the first study of the spectral 
evolution of the continuum and metal line spectrum (see 
Schaerer 1995ab for WR phases). Up to now, theoretical 
studies have either been limited to individual objects, or 
have concentrated on morphological studies of the most 
important P-Cygni lines (Pauldrach et al. 1990). The rea- 
son is directly related to the complex physical conditions 
prevailing in massive early-type stars. This in contrast to 
the situation for intermediate and low mass stars, where 
hydrostatic LTE models provide an good description (but 
see Najarro et al. 1995). To give an example of the status 
of diagnostics in early type stars: at present quantitative 
abundance determinations for OB stars are essentially lim- 
ited to H and He (e.g. Herrero et al. 1992), while the first 
determinations including metals has only been achieved 
recently (Pauldrach et al. 1994). 

The characteristics of our code are described in de- 
tail in Schaerer et al. (1995, hereafter paper I). The code 
takes the non-LTE conditions in a spherically expanding 
atmosphere into account. The photosphere and the wind 
are treated consistently in a unified way, i.e. we do not 
use the core-halo approximation (de Koter et al. 1993). 
In addition we include line blanketing following the opac- 
ity sampling method introduced by Schmutz (1991). A 
description of the sampling technique and a detailed anal- 
ysis of the effects of line blanketing for OB stars is given 
by Schaerer & Schmutz (1994ab, hereafter SS94ab) and 
Schaerer (1995a). 

Our theoretical approach encompasses many aspects of 
the spectral evolution of massive stars. Our model predic- 
tions include the flux distribution from the extreme-UV 
(~ 30 A) to the far-IR (~ 200 /im.) and yield detailed line 
profiles of the strongest H and He lines. We provide ion- 
izing fluxes, which may be of interest for studies of young 
stellar populations, e.g. in clusters, H II regions, and star- 
bursts. Synthetic UBVRIJHKLMN photometric data is 
derived. We present detailed predictions of the UV spec- 
trum taking into account the spectral lines of all elements 
from H to Zn and we compare these spectra with IUE 
observations. 

In Sect. 2 we first present the selected models describ- 
ing the detailed spectral evolution. The continuum spec- 
tral energy distributions, including ionizing fluxes, are dis- 



spectral evolution and compare our results with observa- 
tions. The evolution of the H and He line spectrum, from 
the optical to the IR, is discussed in Sect. 5, where we also 
present an analysis of hydrogen recombination lines which 
can be used as mass loss indicators (Sect. 5.1). The main 
conclusions are summarized in Sect. 6. 
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Fig. 1. a HR-diagram covering the MS phases for initial 
masses of 40, 60, and 85 Mq. The WR stage during the 
H-burning phase of the 85 Mq model is not included. Crosses, 
triangles and stars denote the selected models describing the 
tracks; the square indicates the position of £ Puppis (see Table 
1). b log (/-log T e ff diagram corresponding to the upper panel 



2. Selected models 

In this section we describe the spectral evolution along the 
evolutionary tracks for initial masses M; = 40, 60, and 85 
Mq ■ We have selected several models on each evolutionary 
path, for which we will discuss detailed results. The stellar 
parameters of these models are given in Table 1. Their 
position in the HR-diagram and the log g-log T e ff diagram 
is shown in Fig. 1. In this figure we have included the 
04I(n)f star £ Puppis, analysed in detail by Schaerer & 
Schmutz (1994ab). 

The following entries are given in Table 1: Model num- 
ber (column 1), age (2), present mass (3), the effective 
temperature (4), the luminosity (5), gravity (6), stellar 
radius (7), mass loss rate (8), terminal velocity (9), num- 
ber fraction of hydrogen njj (10), and helium njje (H); 
and the mean molecular weight per free particle \i used to 



Table 1. Summary oi selected models: stellar parameters 
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1 : Schaerer & Schmutz (1994a), model A 



determine the photospheric structure (12). The last col- 
umn gives an approximate spectral classification, which 
has been obtained from a nearest neighbour search in 
the tables of Schmidt-Kaler (1982) using the variables 
(T e ff,Mboi) giving a large weight to the temperature. If 
instead we determine the classification by nearest neigh- 
bour in the HR-diagram, i.e. (log T e ff , log L/Lq), the re- 
sulting differences may reach up to two spectral subtypes 
and two luminosity classes. This should roughly indicate 
the precision of the assigned spectral classification. 

3. Continuous spectrum 

The evolution of the line blanketed flux distribution along 
the MS is shown in Fig. 2 for models along the 60 Mq 
track. The bolometric luminosity increase from the ZAMS 
(model 7) to the TAMS (model 13) is a factor of 1.55. Note 
the progressive appearance of line blanketing features in 
the far UV. This will be discussed in detail in Sect. 4. 

Figure 3 shows a comparison of the EUV to IR spec- 
trum predicted by our non-LTE line blanketed wind mod- 
els (solid line) with a non-LTE H&He plane parallel atmo- 
sphere (Herrero 1994, dotted) and a line blanketed LTE 
model from Kurucz (1991, dashed). The stellar parame- 
ters are chosen as close as possible to that of ( Puppis (see 
Table 1): X , e fj=43 kK, logg=3.7 for the non-LTE plane 



parallel atmosphere, and X , e fj=42 kK, log t/=4.5 for the 
Kurucz model. As expected, the energy distribution of all 
models is virtually identical from the near UV to the near 
IR. At wavelengths A ^ 10 lira, the spherically expanding 
atmosphere shows the characteristic flux enhancement due 
to free-free processes in the wind. The largest differences 
are found shortward of the He II edge (A < 228 A) where, 
compared to LTE and non-LTE plane parallel models, the 
flux is significantly enhanced due to the presence of the 
wind outflow (cf. Gabler et al. 1989). Similarly, non-LTE 
and partly also wind effects increase the flux shortward of 
the He I edge (A < 504 A). Non-LTE and wind effects even 
affect the Lyman continuum flux for lower temperatures. 
These important consequences for the ionizing flux will be 
discussed below. 

3.1. Ionizing spectrum 

In Table 2 we list the predicted number of photons emerg- 
ing at wavelengths shorter than 912, 504, and 228 A re- 
spectively, referred to as qo, q\, and q^- The total ioniz- 
ing luminosity Qi (in photons s _1 ) is obtained by Qi = 
47r(it^i?0) 2 (/j', where R^. is given in Table 1. 

It is difficult to estimate the reliability of the predicted 
ionizing fluxes, since they not only depend on the detailed 
atmospheric modeling (temperature structure, blanketing 
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Fig. 2. Predicted line blanketed flux distributions during the 
MS evolution of a 60 Mq star. Presented are models 7, 9, 10, 
11, and 13 (cf. table 1). The spectra from model 8 and 12 
are very similar to 7 and 13 respectively, and have therefore 
been excluded from this plot. Note that the model close to the 
ZAMS (7) has the largest flux at A < 912 A, and the lowest flux 
longward of the H Lyman edge. With decreasing temperature 
the strength of the blanketing features progressively increases 
in the observable UV spectral range from ~ 1500 to 2000 A 

etc.), but also on the wind parameters (see also Najarro 
et al. 1995). The values for Q2 are affected by the largest 
uncertainty and they strongly depend on the wind den- 
sity. In cases where the radiation shortwards of the He II 
edge originates from deep layers (near to the photosphere), 
the flux at A < 228 A is quite sensitive to the atmo- 
spheric density and temperature structure (cf. also Gabler 
et al. 1992). For example, according to whether the photo- 
spheric structure is calculated taking into account radia- 
tion pressure using OPAL opacities or using only electron 
scattering (cf. Paper I), the predicted He II ionizing flux 
is found to vary by up to 0.6 dex for the models close to 
the ZAMS. This number may be taken as an indication 
of the uncertainty of the Q2 predictions. Another poten- 
tial source of uncertainty for high T e ff models could come 
from the coherent treatment of electron scattering, which 
might modify the He ionization, as recently pointed out 
by Rybicki & Hummer (1994). 

Figure 4 shows the temporal evolution of the ionizing 
luminosity along the tracks. During a large fraction of the 
H-burning lifetime, the ionizing luminosities Qo and Q\ 
remain approximately constant and close to the ZAMS 
value. A constant Qo and Q\ is in fact often assumed in 
calculations of integrated properties of stellar populations 
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Fig. 3. Comparison of emergent spectra from our non-LTE 
line blanketed wind models (solid line) with a non-LTE plane 
parallel model (Herrero 1994, dotted), and a plane parallel LTE 
model of Kurucz (1991, ATLAS 9, dashed line). Also indicated 
are the H Lyman edge, the He I, and the He II Lyman edge, 
as well as the position of some photometric bands (Johnson). 
The stellar parameters are T e jf =42 kK (43 kK for the dashed) 
and log g values of 3.5 to 4.5. The additional parameters for 
the wind model (M, »«,) correspond to £ Puppis (see Table 1) 

(e.g. Vacca 1994). The largest difference is obtained for 
the He ionizing photons along the 40 Mq track, where 
the mean value of Q\ over the MS lifetime is approxi- 
mately half the value at the ZAMS. The situation is sig- 
nificantly different for the He + ionizing photons. As soon 
as the stars have evolved to temperatures T e ff <> 35 kK, 
the wind becomes optically thick to radiation at A < 228 
A, resulting in a rapid drop of the number of He + ionizing 
photons. This corresponds to ages of ~ 2.5 to 3 Myr. In 
a young stellar population the presence of nebular lines 
with ionization potentials above the He + edge thus yields 
a strong upper limit on its age, provided that other ion- 
ization sources such as hot Wolf-Rayet stars or non-stellar 
sources (shock excitation, photoionization by X-rays) can 
be excluded (cf. Schmutz et al. 1992). 

3.1.1. He + ionizing photons 

In Fig. 5 we have plotted the ratio of He II to H ioniz- 
ing photons as a function of effective temperature. The 
results are shown for our line blanketed models. To allow 
for a comparison, we also plot the values obtained with- 
out line blanketing. The spread of the differences between 
blanketed and non-blanketed models is due to differences 



Table 2. Ionizing photon fluxes in cm s from line blan- 
keted models 
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in the calculation of the deep atmospheric structure (see 
above). As demonstrated by Schaerer & Schmutz (1994a) 
line blanketing leads to a higher He ionization and an in- 
creased radiation field at A < 228 A, which explains the 
larger Q2/Q0 ratios obtained for the blanketed models. 
Also given in the figure are the predictions from plane 
parallel non-LTE models for different gravities (Clegg & 
Middlemass 1987). The spherically extended and expand- 
ing atmosphere models predict approximately two to three 
orders of magnitudes more flux in the He II continuum 
than plane parallel atmospheres. As Gabler et al. (1989) 
have shown, this is due to the presence of the wind out- 
flow. Our predictions for Q2 are of the same order as the 
values from the unified models of Gabler et al. (1992). For 
precise comparisons, however, one has to remember that 
for OB stars the He II Lyman flux depends sensitively on 
the atmospheric structure (see above). 

3.1.2. H ionizing photons 

The number of Lyman continuum photons are compared 
with results from previous determinations in Fig. 6. Shown 
is the ionizing photon flux per unit surface area, which 
allows a direct comparison with predictions from plane 
parallel models independently of a radius calibration. We 
compare our models with the widely used results from 
Panagia (1973), the recent calculations from Voels et 
al. (1989) and those from Bohannan et al. (1990). The 
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Fig. 4. Temporal evolution of the ionizing photon fluxes short- 
ward of 912, 504, and 228 A respectively. The short-dashed line 
connects the evolution of the 85 Mq model (crosses), the solid 
line those of the 60 Mq (triangles), and the long-dashed line 
the 40 Mq model (stars). There are three groups of lines for 
Qo , Qi, and Q2 respectively 

latter two results are based on plane parallel non-LTE 
models, which account for backscattered radiation from 
the wind. As expected, for temperatures above T e ff ~ 40 
kK, where the Lyman flux can be well approximated by a 
black-body spectrum, all predictions agree. We add that a 
good agreement is also found with the plane parallel wind 
blanketed models of Kudritzki et al. (1991), calculated for 
35 < T e ff < 51 kK. For the most evolved models, the Ly- 
man continuum flux is, approximately, up to a factor of 
two larger than the Panagia (1973) values — the exact 
value is in particular also dependent on the adopted value 
for M . The enhancement of go is due to the spherical ex- 
tension and non-LTE effects related to the presence of the 
wind. The same result was found recently by Najarro et 
al. (1995), who discuss the physical effect in detail. They 
also point out the possible underestimate of the contribu- 
tion of B stars to the diffuse ionizing radiation. 

3.1.3. He ionizing photons 

Regarding the He ionizing fluxes good agreement is ob- 
tained with the results of Kudritzki et al. (1991, plane 
parallel, wind blanketed non-LTE models): For the range 
covered by their models our predictions lie well between 
their luminosity classes V and I, and differ by less than 
~ 15 %. As mentioned by Vacca (1994), the q\ values ob- 
tained from non-LTE plane parallel atmospheres are sys- 




Fig. 5. Ratio of the He + to H ionizing photons, Q2/Q1, as a 
function of the effective temperature. Circles denote line blan- 
keted models. The effect of line blanketing is illustrated by in- 
cluding the non blanketed models at the same T e jf (same sym- 
bols as in Fig. 4). The square denotes the blanketed £ Puppis 
model from Schaerer & Schmutz (f994a). Solid lines give the 
values obtained from plane parallel non-LTE models of Clegg 
& Middlemass (1987) labeled with their log g value. The large 
increase of the He + ionizing flux is due to the presence of the 
wind 

tematically larger (by factors of 1.6 to 3.7 between T e ff=35 
and 51 kK) than the values derived from Kurucz models 
(see e.g. Fig. 3). For q\, similar words of caution regarding 
the mass loss dependence as for qo are appropriate. 

3.2. Photometric evolution 

We have listed the absolute magnitudes for the Johnson 
UBVRIJKLMN system and the ESO H filter in Table 
3. Following de Koter et al. (1995), we have calculated 
monochromatic absolute magnitudes using the calibration 
of Johnson (1966, F v values from Table IV), Lamia (1982, 
Table 39 for B), and Koornneef (1983) for the ESO filter. 
Note that the calibration of the B filter is quite uncertain, 
which can result in differences up to 0™1 relative to other 
calibrations. To check the assumption of monochromatic 
magnitudes, we convolved our spectra with the response 
functions of Lamia (1982, Table 11, 37 for UBVRIJ). For 
ZAMS models the largest differences were found for RIJ 
(~ 0™ 1-0™ 15 smaller). Differences in U can amount to 
0™1, while B and V show smaller differences. Therefore, 
for more precise determinations the appropriate response 
functions should be used. 



Fig. 6. Logarithm of the Lyman continuum flux (photons 
cm -2 s -1 ) versus effective temperature. The predictions from 
our line blanketed atmosphere models are indicated as crosses 
for the 85 Mq track, triangles (60 Mgtrack), and stars (40 
Mq track). The solid lines shows the relations from Panagia 
(1973) for his ZAMS and luminosity class I. The dashed line 
connects the wind blanketed models from Voels et al. (1989) 
and Bohannan et al. (1990) (circles). 

In particular our predictions of infra red (IR) magni- 
tudes, taking the effect of the spherically expanding non- 
LTE atmosphere into account, may provide a useful base 
for many applications in studies of hot star populations. 
However, one should keep in mind that the IR flux de- 
pends on the adopted wind parameters, in particular on 
the wind density. As an example: for the supergiant model 
12 an increase of M jv^ by a factor of 3.5 (8) yields fluxes, 
which differ by less than 0™05 (~ 0™1) for UBVRI, by 0™1 
(0™5) for K, and by up to 0™6 (1™8) for the N filter. Our 
predictions should therefore be used with caution for stars 
with significantly different mass loss rates and/or wind ve- 
locities. For LBV's the influence of parameter changes on 
the emergent spectrum and on photometric quantities has 
been studied by de Koter et al. (1995). 

3.3. UV continuum indices 

The optical continuum of hot stars carries very little infor- 
mation about gravity and temperature because it is in the 
Rayleigh-Jeans domain of the spectrum. For photometric 
studies one is therefore led to use bands located in the UV 
spectral range. 

In Table 4 we list the color indices from the contin- 
uum flux of our line blanketed models for UV passbands 



Table 3. Monochromatic absolute magnitudes lor the Johnson UBVR1JKLMN system. Also given is the LSO H magnitude 
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Fig. 7. UV color-magnitude diagram for the (1400 — V) index. 
The absolute visual magnitude is taken from Table 3. Using 
UV colors, different evolutionary phases are clearly separated 
(symbols as in Fig. 1). The dashed and the dotted line show 
the relations for the mean values for dwarfs and supergiants 
respectively given by Fanelli et al. (1992). The large stretch 
in (1400 — V) illustrates the power of UV photometry as a 
diagnostic tool for hot stars 

as defined by Fanelli et al. (1992, hereafter FCBW). The 
UV bands are rectangular boxes of width AA = 200 A, 
approximately centered at the given wavelength (in A). 
Using IUE low resolution spectra FCBW have built mean 
spectra for a large number of groups of spectral types and 



luminosity classes (LC), which were then used to derive 
UV colors and line indices. We will compare our calcula- 
tions with their results. 

In Fig. 7 we plot our MS tracks in the M v vs (1400-U) 
diagram. The (1400 — V) color index varies by up to two 
orders of magnitude between spectral types ~ 05 V and 
B2 I. This clearly illustrates the usefulness of UV indices 
as a diagnostic tool for hot stars. The mean observed val- 
ues from FCBW for dwarfs and supergiants, up to the 
spectral type 04, are also shown in Fig. 7. Our models 
cover the entire observed range between LC V and I, the 
most extreme models being more luminous than the mean 
LC I. In Fig. 7 our ZAMS models appear to be ~ m l-0 m 2 
too blue in (1400 — V). However, the standard deviation 
(using equal weights) from the stars in the 03-6 V group 
of FCBW (see Fig. 7) is m 23. Thus, our results for the ob- 
served continuum slope between 1400 and 5500 A agrees 
within one sigma. We add that our predicted (1400 — V) 
colors also agree with those derived from Kurucz (1991) 
models. 

To make a more detailed comparison of the pre- 
dicted UV continua with observations we have plotted 
the (1400 - 1700) vs. (1700 - 2600) color-color in Fig. 8a. 
Again, the mean color indices from FCBW for dwarfs and 
supergiants of spectral types between 04 and B9 are su- 
perposed in this figure. As the majority of our models 
are roughly of spectral types 03 to 09, they are crowded 

on the left of the figure, where (1400 - 1700) m 5 and 

(1700-2600) l m 4 to -l m 5. This is in good agreement 

with the observations. Note also that for these models 
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Fig. 8. a UV color-color diagram for the bands from FCBW defined in the text. Our models are shown by the same symbols 
as in Fig. f . Note that the models of approximate luminosity class V, and IV (cf. Table I) are indistinguishable and located at 

(1400-1700) m 5 and (1700-2600) f m 4 to -f m 5. The giant and supergiant models are well separated. The long-dashed 

line indicates the position of our models if we exclude line blanketing. Observations (mean colors from Fanelli et al. 1992) are 
shown for Dwarfs (dashed) and Supergiants (dotted). The labels along these mean relations indicate the position of the spectral 
type groups, which have been used to derive the relations. The most evolved models are found to have an excess of up to ~ 
m I5 mag in the 1700 band. 

b Same as left panel showing (1400 — 2600) vs. (2600 — V). The models are in good agreement with the mean observed relation. 



the predicted blanketing effects lead only to small changes 
of the colors considered above. Surprisingly however, the 
remaining models, corresponding to later types and su- 
pergiants, are located outside the region spanned by the 
observations. The systematic shifts indicates mainly that 
the observed (1400 — 1700) color must be bluer than the 
models predict. A comparison with the (1400 — 2600) vs. 
(2600 — V) diagram (Fig. 8b) shows a good agreement with 
the observations. This implies that the difference noted 
above must essentially be due to the 1700 band. Indeed, 
for the most discrepant supergiant model (model 6) a de- 
crease of the flux in the 1700 band by only 0™ 15 suffices 
to reconcile our models with the observed (1400 — 1700) 
and (1700 — 2600) colors for supergiants. 

To investigate the cause of this difference, we compared 
our models with several IUE spectra of supergiants, as 
used by FCBW. This allows us to explain the discrepancy: 
The observed spectrum plotted in Fig. 9 (HD 165024, B2 
lb) is from the BO-2 I group. Contrary to our model, the 
observed spectrum has a steeper gradient between the 
1400 and 1700 bands. As explained above this difference 
is essentially due to a larger flux depression in the 1700 
band. The most likely explanation is that over a large frac- 
tion of the 1700 band numerous Fe IV lines (cf. Bruhweiler 



et al. 1981) provide important blocking, which is underes- 
timated in our models. In fact, in this wavelength domain 
the supergiant spectra look very similar to the spectra of 
Of/WN objects, for which Schaerer (1995a) demonstrated 
that the large flux depression is mainly due to Fe IV. A 
more quantitative discussion of several line blanketing fea- 
tures is given below. We thus conclude that in the 1700 
band our models yield an overestimate of up to 15 % in the 
flux, while for the other far-UV colors we find a general 
agreement with the observations. 

4. Evolution of the UV line blanketed spectrum 

Our line blanketed atmosphere models allow us to pre- 
dict a large number of observable line features, which are 
very well suited for comparisons with IUE low resolution 
spectra. Before doing so, we again (see Paper I) like to 
stress that three important assumptions enter in the cal- 
culation of spectra by Monte-Carlo simulation. First and 
second, those on the ionization and excitation structure, 
as discussed in SS94ab. Third, those on the line forming 
mechanism. In the present work we concentrate on wind 
features since photospheric broadening is not treated. In 
our Monte-Carlo method, we include line scattering in the 
Sobolev approximation taking the continuum opacity into 



Table 4. UV continuum color indices (in magnitudes) 
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account. To be precise, collisional deexcitation and the ef- 
fect of branching are not taken into account. The weakest 
point in our calculations is likely the neglect of line broad- 
ening yielding only a poor treatment of photospheric lines. 
We would like to stress that the problem of non-LTE line 
blanketing in hot star atmospheres is a very complicated 
one and a full treatment is as yet not feasible (see Paul- 
drach et al. 1994). Our calculations are a first step towards 
this goal. 

Due to the adopted frequency binning, the spectral 
resolution varies nonuniformly. It is largest longward of 
the H Lyman edge at 912 A and the He II Paschen edge 
at 2046 A. The resolution ranges between A/AA ~ 100 - 
1000 with a typical value of 500. At all wavelengths the 
relative statistical error of the predicted flux is less than 
10" 3 . 

In the following, we will study the evolution of the most 
important wind lines, but we will also discuss other metal 
line features. The predicted UV spectral evolution for the 
Mi = 40, 60, and 85 M Q tracks is shown in Figs. 10, 11 
and 12 respectively. The illustrated spectral range reaches 
from ~ 900 to 2200 A. The main sequence evolution ap- 
proximately starts with spectral types 05 to 03 dwarfs, 
which then evolve to B2 to 05 supergiants (cf. Table 1). 
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1500 2000 2500 

A [A] 

Fig. 9. Comparison of the predicted UV spectrum of model 6 
(solid line) with IUE low resolution observations of HD 165024 
(B2 lb, dashed). The observed spectrum has been dereddened 
using E(5-U)=0 m l (FCBW) and the Seaton (1979) law. The 
scaling of both observations and the model spectrum are arbi- 
trary. Also shown are the positions of the broad far-UV bands 
1400, 1700, and 2600. The thin solid line gives the slope of the 
model spectrum between these bands (drawn by hand). The 
thin dashed line is obtained by shifting the (1400 — 1700) and 
(1700 — 2600) slopes of the model to the observed spectrum. 
The observations clearly show a bluer (1400 — 1700) index prob- 
ably due to important blending in the 1700 band 



CoStar: 40 ML, Z = Q.Q2, main sequence 
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Fig. 10. Synthetic UV spectra showing the spectral evolution on the 40 Mq track (models 1-6). Plotted is the logarithm of 
the emergent luminosity. The evolution along the MS proceeds from top (dwarf) to bottom (supergiant). Starting with the 
second model, each spectrum has been shifted downwards by 0.7 dex with respect to the previous one, in order to allow a good 
comparison. The marks on the top and bottom indicate the location of the CNO and Si lines taken from the lists of Bruhweiler 
et al. (1981) and Dean & Bruhweiler (1985). The strongest CNO, and Si features are labeled. Note also the strong Fe III 1920 
A feature in supergiants 




Fig. 11. Same as Fig. 10 for the 60 M Q track (models 7 to 13) 




Fig. 12. Same as Fig. 10 for the 85 M Q track (models 14 to 16) 



Ovi AA 1032, 1038, N v AA 1239, 1243, Si iv AA 1394, 
1403 and C IV AA 1548, 1551 are clearly present in the 
synthetic spectra plotted in Figs. 10 to 12. A first inspec- 
tion of the strong wind lines yields the following: 

(a) The predictions for the N V resonance line follows 
the observed decrease in line strength towards later spec- 
tral types. The spectra of the B type stars, show some 
tendency to produce less N V than observed, which re- 
flects the problem of super-ionization. The observed N IV 
A 1718 line shows no clear trend in the domain of our 
calculations (Snow et al. 1994). The observed profiles are 
of different types including photospheric lines, shifted or 
asymmetric absorption and P-Cygni profiles. We find the 
same variety of profiles in our models. 

(b) The strong luminosity dependence of the Si IV line 
is reproduced. 

(c) The observed C IV line in 03-5 dwarfs is very 
strong. Our ZAMS models, however, do not produce this 
feature. We will return to this problem later on. For the 
domain covered by our models, the observations of the 
C III A 1175 line show a gradual transitions from pho- 
tospheric to P-Cygni profile (Snow et al. 1994). The ap- 
pearance of this C III feature in our models follows the 
observed trend. 

(d) We predict a rather strong O V A 1371 P-Cygni 
line for the hottest models, while it is mostly observed in 
absorption. This might be due to an overestimate of either 
the ionization fraction of O V, or the excitation of the cor- 
responding level. We note that for the hottest models, O V 
is indeed predicted to be the dominant ionization stage in 
contrast to the calculations of MacFarlane et al. (1993) 
and Pauldrach et al. (1994). Similarly the overionization 
of O is also apparent if we compare the O IV A 1338 and 
O V A 1371 lines with the observations compiled by Snow 
et al. (1994). 

Let us now make a systematic and quantitative com- 
parison of the predicted spectra with observations. We 
again use the results of FCBW, who introduced several 
far-UV spectral line indices which have been derived for 
a large number of stars using IUE low resolution spectra 
covering luminosity classes V to I and spectral types 03 
to K or M. 

In Fig. 13 we illustrate the position of the narrow 
bands of FCBW (bands 2-20), and of three bands de- 
fined by Nandy et al. (1981, called 54, 56 and 57). These 
bandpasses were chosen to include the most prominent 
UV absorption and P-Cygni features and adjacent "con- 
tinuum" regions. The characteristics of the bands 54, 56, 
and 57 are given in Table 5. Spectral indices, derived from 
combinations of these bandpasses, quantify the strength of 
spectral lines. Line indices are determined by mean fluxes 
over three bandpasses, one being centered on the wave- 
length of the feature, and two on bracketing sidebands. 
All indices are measured in magnitudes. The indices of 
FCBW are calculated according to J; = -2.blog(F x / F x ), 
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Fig. 13. Detail of normalised UV spectrum (model 10) to il- 
lustrate the position of the far U V narrow bands listed in Table 
5. Bands 2-20 are from FCBW; 54, 56 and 57 from Nandy et 
al. . The position of two of the broad bands of FCBW is shown 
at the bottom. As in Fig. 10 the marks indicate the positions 
of observed CNO and Si lines 



where F x is the mean flux in the feature band i, and F x is 
the linearly interpolated continuum flux expected at the 
line center, computed from the fluxes in the bracketing 
sidebands. An index I > thus indicates a net absorp- 
tion. Table 6 lists the line indices, defined by their central 
bandpass and sidebands. Note that the C IV P-Cygni line 
index is defined by I = — 2.5 log(.Fio/-fi2); where Fio and 
F12 are the mean fluxes in the bands 10 and 12 respec- 
tively. 



Table 5. Characteristics of the far-UV bands defined by Nandy 
et al. (1981). The first column gives the identifier of the band 
(cf. Fig. 13). The initial and final wavelength, and the total 
width of the bands (in A) are given in cols. 2, 3, and 4 respec- 
tively. Column 5 indicates whether the band is a continuum 
sideband (C) or the actual bandpass with the line feature (L). 
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Following Nandy et al. (1981), the S1920 line index is 
determined by k = -2.5 \og{2F[/ (F[ + F{*)), where F l x 
and F x are the total fluxes in the left and the right side- 
band respectively (cf. Table 5). The UV line indices de- 
rived from our models are plotted in Fig. 14 as a function 
of the color (2600 — V), which serves as a temperature 



Table 6. iar-UV spectral line indices 



Feature name Central wavelength Bandpass Sidebands 
[A] 



Indices from Fanelli et al. (1992): 
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indicator and allows a direct comparison with the obser- 
vations given by FCBW. 

4.1. Quantitative comparison of CNO and Si lines 

First we shall study the four indices which describe the 
C IV resonance line. The first three indices (centered at 
1540, 1550, and 1560 A) are measuring the strength of 
the blueshifted absorption, the central absorption, and the 
redshifted emission respectively. The fourth or 'P-Cygni' 
index measures the ratio of the emission to the absorp- 
tion, characterizing the net strength of the wind feature. 
The most sensitive indices to quantify the P-Cygni line are 
clearly the blueshifted absorption and the P-Cygni index, 
which show variations of up to 1 mag over a small range 
of (2600 — V). As is clear from Fig. 14, our models re- 
produce both the magnitude and the trends of these C IV 
indices well. However, as already mentioned above, the 
ZAMS models predict a too weak C IV A 1550 line. This 
is clearly an indication that in the hottest models, where 
C V is predicted as the dominant ionization stage, the car- 
bon ionization seems to be overestimated. For the coolest 
supergiant models we predict slightly less blueshifted ab- 
sorption (1540) and too much redshifted emission (1560). 
This may also be due to an overestimated carbon ioniza- 
tion balance. 

For silicon, FCBW introduced the Si IV 1397 index po- 
sitioned to measure the absorption trough of the Si IV res- 
onance doublet. Its strong luminosity dependence is well 
known (e.g. Walborn k Panek 1984). Figure 14 shows a 
systematic shift of the Si IV index. The models have sys- 
tematically too much flux in the Si IV band with respect 
to the sidebands. This has to be taken with caution since, 



diet a too strong O V A 1371 line, which reduces the flux 
in sideband 4 used to derive the Si IV 1397 index. Artifi- 
cially eliminating the O V blueshifted absorption leads to 
an increase of the Si IV 1397 index by ~ 0™15-0™2 for the 
hottest models. In addition, one must consider that for 
the early type main sequence stars which do not show the 
Si P-Cygni feature, the 1397 band is also strongly pop- 
ulated by numerous Fe V blends (Bruhweiler et al. 1981). 
This may also explain part of the shift between our models 
and the observations. In view of the above, we conclude 
that the shift in the Si IV 1397 index is most likely not re- 
lated to an error in the prediction of the silicon ionization 
structure. 

A comparison of the actual profile shape of Si IV with 
the atlas of Snow et al. (1994) indeed reveals that the ob- 
servations are well reproduced. Especially the tempera- 
ture, below which the Si IV line shows a P-Cygni profile 
is correctly predicted for each track. Moreover, the trend 
of an increased strength of the absorption trough of the P 
Cygni line with decreasing temperature is recovered. 

The BL 1302 feature is produced by several lines of 
Si II, Si III, and O I (Fanelli et al. 1992, Shore k Sanduleak 
1984), while BL 1425 is due to a blend of Fe V, C III, and 
Si III. Given the fact that several elements contribute to 
the corresponding blends, there seems no simple interpre- 
tation to explain their behaviour and the small differences 
with our models. The remaining indices (BL 1664, 1719, 
and 1853) show no important trends for the hot stars con- 
sidered in this work. The model predictions for these in- 
dices are in tentative agreement with the observations. 

4-2. Quantitative comparison of Fe features 

Of particular interest are line indices which are indica- 
tors of the strength of iron lines, since these lines are the 
most important contributors to blanketing as well as to 
the radiative force (e.g. Abbott 1982). The models show a 
strong increase of the metal lines along the evolution from 
the ZAMS towards the cooler supergiant phase. A very 
prominent Fe III feature appears at A ~ 1920 A in the late 
O to early B supergiants. This feature will be investigated 
in detail in Sect. 4.2. First, we discuss the diagnostics of 
Fe iv k v. 

The Fe V 1453 index shown in Fig. 14 is a measure 
of a blend of several Fe V lines (cf. Nandy 1977). Its ob- 
served strength is 0™2. As shown by the observations 
of FCBW, and apparent in Fig. 14, this index shows a 
strong temperature dependence. While our coolest mod- 
els approximately agree with the observations, we seem 
to predict too weak a Fe V 1453 absorption for the less 
evolved models. However, as for the Si IV line discussed 
above, this index is very sensitive to the predicted O V 
line strength: By artificially eliminating the O V contribu- 
tion we find that the indices which use sideband 4 (among 
which Fe V 1453) may be too low by ~ 0™1. Correcting 
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Fig. 14. UV line indices in magnitudes as a function of the (2600 — V) color, which serves as a temperature indicator. The 
symbols (cf. Fig. 1) denote the values predicted by our models. The dashed and the dotted line indicate the mean values derived 
by FCBW for Dwarfs (LC V) and Supergiants (I) respectively 



for this effect a difference of up to 0™1 remains. We there- 
fore conclude that our Fe V 1453 index may point to a 
somewhat underestimated Fe V ionization. 

Another of the FCBW indices measures the contribu- 
tion of several Fe lines. The BL 1617 absorption is at- 
tributed to a blend of Fe IV and Fe V lines (cf. Dean & 
Bruhweiler 1985). The observations show a clear luminos- 
ity separation. Our models predict a behaviour similar to 
the Fe V 1453 index. 

4.2.1. The Fe III feature at 1920 A 

As mentioned above, a striking metal line feature around 
1920 A appears in our late O and early B type models. 
This broad absorption feature, which we shall refer to as 
the '1920 feature', is due to important blocking of several 
Fe III lines. It has already been identified in TD1 observa- 
tion of B type stars (Thomson et al. 1974, Nandy 1977). 
Subsequently the 1920 feature was shown to increase from 
late to early B stars (Swings et al. 1976), and to decrease 
again for O stars (Nandy et al. 1981). To quantify this 
absorption feature we have computed a photometric mea- 
sure of its strength, S1920 defined by Eq. 4 and Table 5. It 
measures the depression in a 30 A band centered at 1920 
A with respect to continuum points chosen at 1700 and 
2200 A. 

In Fig. 15a, we plot the S1920 index derived from our 
models as a function of effective temperature. In our mod- 
els the 1920 feature appears first for temperatures T e ff ;S 



33 kK. For these temperatures, we find a strong correlation 
of its strength with T e ff . Down to our coolest model with 
T e ff ~ 18 kK, the S1920 index increases steadily by up 
to 0v Q 6-0v Q 8 with decreasing temperature, as expected by 
the progressive decrease of the Fe ionization. Note, that 
in Fig. 15a we have added the value of the S1920 index 
derived from a model of the Galactic Center He I emis- 
sion line star AF (Ofpe/WN9; parameters from Najarro 
et al. 1994, model Bif). It shows the strongest 1920 feature 
of all models calculated so far, and appears as a natural 
continuation of the early B hypergiants. First, this indi- 
cates that if the derived stellar parameters are correct, the 
effects of line blanketing are expected to be particularly 
strong. Secondly, this may also be an additional hint show- 
ing the intimate relation between LBV's (cf. P Cygni) and 
Ofpe/WN9 stars. Both points will be discussed in a forth- 
coming paper. 

To allow for a comparison between the predicted Fe III 
feature and observations we have plotted the S1920 index 
as a function of absolute visual magnitude in Fig. 15b. 
The corresponding observations are given by Nandy et 
al. (1981, their Fig. 8). We note that the strength of the 
predicted 1920 feature agrees well with the observations, 
and we correctly reproduce the observed trend of increas- 
ing strength from O to early B stars. 

To compare the predicted T e ff trend with observa- 
tions, we have derived the S1920 index from IUE spec- 
tra for several extensively studied dwarfs and supergiants 
with relatively low reddening and well determined tern- 



o 

05 



0.2 




0.8 



0.6 - 



0.4 



0.2 



0.E 



o 

05 



0.4 



0.2 



4.6 4.7 



0.2 



1 


1 1 1 1 1 


1 1 1 1 | 1 1 1 1 


1 


- 






b 


- 


# 




- 




A 




- 


- 


A 




- 


- 




* 


- 












A 








Tr 








A 




_ 




X X 


_ 






A 








X 




- 

- 1 — 


I...! 


A 

i i i i 1 , , , 


- 




a 
O 


1 D 


u 






M v 





Fig. 15. a Predicted line index (in mag) measuring the Fe III 1920 A feature as a function of logT e jf. The same symbols as in 
Fig. 1 are used. Also shown is one model for the Of/WN Galactic Center star (AF). The dashed line connects observations for 
stars with well determined T e jf . Note the strong temperature dependence of S1920, which is confirmed by the observations, 
b Same as panel a, but shown as a function of Mv- This figure can directly be compared with observations of OB supergiants 
of Nandy et al. (1981) 



peratures. All spectra were dereddened with the Seaton 
(1979) law using the following color excesses. For 9 Sgr: 
E(5 - V) = m 33 (Fanelli et al. 1992), ( Ori A: E(5 - V) 
= m 06, k Ori: E(B-V) = m 05 (both Bastiaansen 1992), 
and P Cyg: E(5 - V) = m 43 (Lamers et al. 1983). The 
temperatures are from Bohannan et al. (1990), Voels et 
al. (1989), Lennon et al. (1991), and Lamers et al. (1983; 
cf. Pauldrach k Puis 1990). Since it might be argued that 
the luminous blue variable P-Cygni is a special case, which 
may not be well represented by our models, we have ex- 
tended our comparison at low T e ff by adding observations 
of objects that have no individual T e ff determination. We 
use two B2Ia supergiants \ 2 Ori an(1 10 Per. From the 
spectral type, Lennon et al. (1993) assign T e ff ~ 17.5 kK 
to these stars, while the Schmidt-Kaler (1982) calibration 
gives T e ff = 18.5 kK. The E(5 — V) color excesses were 
obtained from Geneva photometry (Rufener 1988) and the 
intrinsic colors from Cramer (1993). The resulting S1920 
indices are 0™49 and 0™50 mag for \ 2 Ori an d 10 Per re- 
spectively, which agrees well with the value determined for 
P-Cygni. The observed S1920 indices, plotted in Fig. 15, 
are found to confirm the predicted strength and the tem- 
perature dependence of the Fe III feature. 

Apart from showing a temperature trend, one could 
also expect the Fe III blanketing feature to depend on the 
wind density. We note however, that although the models 
with the strongest S1920 index have terminal velocities 



roughly a factor of two larger than those of \ 2 Ori and 10 
Per (Prinja et al. 1990), this effect is compensated by the 
larger mass loss rate of the models. The comparison of the 
predicted 1920 feature with the observations is thus well 
posed. 

The above results show that the calculated ionization 
of Fe agrees well with the observations. Since Fe is the 
dominant contributor to line blanketing this provides a 
strong support for our models. In conclusion, the 1920 A 
feature provides an excellent test for atmosphere models 
at T e ff 33 kK. For hotter stars however, where higher 
ionization stages of Fe are dominant, we must conclude 
that it is more difficult to obtain constraints of similar 
quality on the behaviour of Fe. 

5. Evolution of the optical and IR line spectrum 

The most important spectral signatures of OB stars in 
the optical and IR are the hydrogen and helium lines (see 
Walborn k Fitzpatrick 1990, Hanson k Conti 1994). To 
illustrate the predictions made by our models we present 
in Fig. 16 the evolution of some important H and He lines 
in the 3500 A to 4.4 //m interval for the H-burning phase 
of the 60 Mq model. The lines have been obtained from 
the line blanketed non-LTE calculations assuming a total 
thermal plus turbulent broadening of 15 km s _1 . In the 
present treatment, our predictions should be reliable for 
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Fig. 16. a and b. Predicted line blanketed non-LTE spectral evolution of some important hydrogen and helium lines of a 60 
Mq star during the MS. Lineprofiles have been superposed to the continuum fluxes. The units and the models shown are the 
same as in Fig. 2. a Visible spectrum. The strongest line plotted is Ha (A 6562.8 A), which changes from absorption to emission 
with increasing mass loss going from the ZAMS to the TAMS. 

b IR spectrum covering the IJKL bands. The strongest lines are Paa A 1.8751 /tin, Br/? A 2.6252 /tin, and Bra A 4.0512 /tin, 
which show a similar behaviour as Ha depending on the mass loss rate 



wind lines, while pure absorption lines are less reliable 
because of the use to the Sobolev approximation and the 
neglect of pressure broadening (cf. de Koter et al. 1993). 
Lines contributing to blends are plotted separately. The 
absence of certain lines is due to the choice of the atomic 
models. 

The most striking change in the optical spectrum is the 
strong increase of the Ha emission towards the end of the 
MS evolution. This basically reflects the increasing mass 
loss (see Paper I), which conversely can be used to deter- 
mine mass loss rates from Ha (see below). The strongest 
lines in the IR region shown here are from the Paschen and 
Brackett series, which show a similar increase with mass 
loss as Ha. This will be discussed in Sect. 5.1. Essentially 
all predicted Bra/Paa equivalent width ratios (neglecting 
blends; cf. below) are between 2.8 and 3.8. Only for the 
models with the weakest emission this ratio is found to 
increase considerably (up to 8). 

Of particular interest is also the spectrum in the K 
band. In this spectral region hot stars in highly obscured 
star forming regions can be observed. A preliminary spec- 
tral classification of OB stars in the K band has been 
presented by Hanson & Conti (1994). Their observations 
reveals four main atomic transitions, due to hydrogen, he- 
lium, and nitrogen. Here we present the predicted H and 
He lines in the 2 to 2.25 jj.ni interval, i.e. Br7 (calculated 



without the He I (7-4) blend), He I 2.058, and the (10- 
7) He II 2.1885 transition. The observed He I triplet line 
at 2.113 was not included in our line blanketed non-LTE 
calculations. The theoretical line profiles obtained for the 
main sequence evolution of the 40 and 60 Mq track are 
plotted in Fig. 17. Since pressure and rotational broaden- 
ing are not taken into account in the present calculations 
the absorption lines are only indicative. We simply plot a 
triangle profile with a constant width (±0.1 in At; /vo) and 
the relative flux at line center for pure absorption lines. 
Their interest is to illustrate the relative changes during 
the main sequence evolution. 

Qualitatively, the evolution of the illustrated lines is 
identical on the two tracks. With the progressive temper- 
ature decrease and the mass loss increase, the absorption 
in the He I and the Br7 line gets progressively weaker and 
switches to emission for the most evolved MS models. The 
He II line is very weak. It is strongest in ZAMS models, 
where it is found in absorption. 

A comparison with the observations from Hanson & 
Conti (1994) yields the following: The strongest line in 
dwarf models (1 & 2 and 7 & 9) is Br7, He II 2.189 may 
show absorption while He I 2.058 is very weak or absent. 
This qualitatively agrees with the observations of lumi- 
nosity class V objects. Our supergiant models (5 & 6 and 
life 13) show for the Br7 and the He I line a net emis- 



TJ 2 
0) 



_ 4 



_ 3 



Hel 2.058 fj. 



'I AMS 



ZAMS 



Br 7 IIcII 2.189 



2.05 2.1 2.15 2.2 



2.25 



T3 2 
0) 



Hel 2.058 



1 3 



_ 11 



_ 10 



TAMS 



ZAMS 

i 



Br 7 Hell 2.189 



2.05 2.1 2.15 2.2 

A | fjum | 



2.25 



Fig. 17. a Predicted K band line profiles for the main sequence evolution of the 40 and b 60 Mq tracks. Shown are the He I 
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The results for model 8 and 12, which are not shown, are very similar to 7 and 13 respectively 



sion with an equivalent width ratio of Br7/He l~ 0.5 to 4. 
The strongest emission is W eq ~ -13. and -4.4 A for Br7 
and He I 2.058 respectively. From these emission lines the 
model predictions close to the TAMS are thus more simi- 
lar to late Olf types than to early B supergiants (see Han- 
son & Conti). Although very weak, the He II 2.189 line is 
predicted in emission contrary to the observations show- 
ing He II in absorption in Of supergiant spectra. We have 
assigned spectral types of ~ BO-2 I to the TAMS stars, 
using a nearest neighbour search in the T e ff-My M \ plane. 
However, if one performs a search in logT e ff- log_L, one 
obtains spectral type 071 to 09.51. The comparison with 
late Olf stars is therefore justified. 

We conclude that the discussed K band lines of H 
and He are in tentative agreement with the general ob- 
served behaviour. This illustrates the capability to repro- 
duce near IR lines with our non-LTE models. More de- 
tailed work will be required to allow determinations of 
stellar parameters of OB stars, based on IR spectra only. 

5.1. Hydrogen recombination lines as mass loss indicators 
5.1.1. Optical spectrum: Ha 

The mass loss rates of hot stars may be determined using 
the Ha recombination line. Given the strong Ha emission 
in hot stars, this technique can readily be applied to extra- 
galactic objects that are not detectable with other tech- 
niques such as radio or sub-millimeter measurements. The 
method has been developed by Leitherer (1988), Scuderi 



et al. (1992), and Lamers & Leitherer (1993, henceforth 
LL93). 

The procedure can be summarised as follows: Provided 
the Ha luminosity _L(Ha), i.e. the distance is known, one 
subtracts a certain photospheric contribution from the ob- 
served Ha. Assuming an optically thin line and an isother- 
mal wind one can then derive the mass loss rate pro- 
vided additional stellar parameters (i.e. radius, temper- 
ature, wind velocity) are known. An uncertainty in de- 
riving _L(Ha) is the assumption of an underlying photo- 
spheric absorption profile. The unknown strength of the 
He II A 6560 blend (cf. LL93), which also contributes to 
the final observed Ha line may also introduce a source 
of error. More refined procedures, which avoid these as- 
sumptions were recently presented by Puis et al. (1995). 
The important drawback of all these methods is that they 
require knowledge of photospheric parameters (in any case 
T e flf, \ogg), distance (to derive R^) and terminal velocity 
(e.g. from UV observations). 

For applications where the above requirements cannot 
be fulfilled it may be useful to provide a simpler procedure 
to derive M . To this end, we first present a straightforward 
fit formula to our atmospheric model results, which allows 
to derive the mass loss rate from a sole Ha emission equiv- 
alent width measurement. Such a relation is particularly 
useful for stars with unknown distance and/or unknown 
wind velocities. In a second step we present an improved 
relation, taking into account the stellar parameters. 
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Fig. 18. a Mass loss rate as a function of predicted log IT(Ha) (emission) including the He II A 6560 blend (same symbols as 
in Fig. 1; 3 additional models, cf. text). The solid line is a fit to the data (Eq. 1), the dashed line is the relation from Gabler et 
al. (1990). 

b Generalised curve of growth for Ha obtained from our models (same symbols as in panel a). Shown is Q 2 as a function of 
log IT(Ha) (emission). The solid line shows a linear fit to the data, which is used to derive the coefficients a and b used in Eq. 3. 
Squares show the results from the unified models of Puis et al. (1995, Table 3), which are fitted by the dashed line 



In Fig. 18a, we have plotted the mass loss rate as a 
function of the predicted equivalent width of Ha. In ad- 
dition to the models from Table 1 we have included three 
extra models (found at 0.2 < log IT(Ha) < 0.4 A) to ob- 
tain a more regular distribution of IT(Ha). Since our mod- 
els treat the entire atmosphere, including the photosphere 
and the wind, the predicted Ha flux can directly be com- 
pared to the observations without having to recur to any 
arbitrary assumptions about a photospheric contribution 
and the He II blend. Note that we have only shown the 
predictions for models having net Ha emission, since for 
weaker lines the present treatment may underestimates 
the photospheric absorption, and is therefore not suffi- 
ciently reliable (cf. de Koter et al. 1993). A least square 
fit to our data yields: 

logM = (0.568 ± 0.100) loglT(Ha) - 5.647± 0.093 (1) 

with a residual rms of ~ 0.17 dex. Note that we obtain a 
quite similar behaviour as given by the relation of Gabler 
et al. (1990, Eq. 4), which is also plotted in Fig. 18. 

As is immediately apparent in Fig. 18a, IT(Ha) is 
also dependent on other basic stellar parameters than 
mass loss. For given M the Ha emission is stronger in 
less massive stars. Given their lower luminosity the same 
M is found at lower T e ff, which explains the major part 
of the increased Ha emission. Puis et al. (1995) found 



that the Ha equivalent width is expected to scale with 
Q 2 = M 2 R~ 3 Vcx^^ 3 ct, where ct is a 'temperature correc- 
tion factor' with respect to a given reference temperature. 
Following these authors we approximate ct by: 
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where cl = 5.184, c2 = 2.337, and T e ff is expressed in 10 4 
K. In Fig. 18b, we plot the behaviour of Q 2 as a function 
of IT(Ha). This shows that at least half of the scatter 
present on the left panel can be explained by the scaling 
relation for optically thin emission taking the stellar pa- 
rameters R^, T e ff and into account 1 . The best fit of 
logQ 2 vs. loglT(Ha), given as the solid line in Fig. 18b, 
can be inverted to derive M . This results in: 
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log W + i log (Rlv^cr 1 ) 
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where W is in A, R* in units of Rq, v^, in km s 1 , and 
ct is given by Eq. 2. This equation defines the rescaled 



1 The apparent scatter in our data stems from models with 
Teff^ 32 kK, which have considerably larger 63 departure co- 
efficients accounting for their 'excess' emission. 
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Fig. 19. a Mass loss rate as a function of the rescaled Ha equivalent width W' for our models (symbols as in Fig. 18a) The 
solid line shows the relation given by Eq. 3. 

b Comparison of mass loss rates (Mat) using Eq. 1 with M determined by Puis et al. (1995, Table 11, col. 3). Plotted are all 
stars compiled by Puis et al. showing net Ha emission (see LL93). The thin dotted line is the one-to-one relation. Stars show 
the individual values obtained from Eq. 1 (mean relation: solid line). Squares denote the values obtained from Eq. 3 (mean 
relation: dashed line) 



equivalent width W' which also depends on the coefficient 
a. For Ha we obtain a = 0.923 and b = —21.122. Equation 
3 thus allows a more accurate derivation of M from the 
Ha equivalent width provided R^, T e ff and are known. 
The corresponding relation between W' and M is shown 
in Fig. 19a. 

In Fig. 18b, we have also plotted the results from the 
unified models from Puis et al. (squares). A fit to their 
data (dashed line) yields a = 0.704 and b = —20.737. Sur- 
prisingly, the comparison with our models shows a sys- 
tematic difference, the models from Puis et al. predicting 
lower equivalent widths for identical stellar parameters 
(i.e. same ordinate Q 2 ). Line blanketing, which is the main 
additional ingredient in our models with respect to Puis et 
al. only accounts for small changes. This difference must 
be due to (a) differing atmospheric structures in the tran- 
sition zone between photosphere and wind, which sensi- 
tively affects the Ha equivalent width 2 (cf. SS94a), or (b) 
the neglect of pressure broadening, which underestimates 
the remaining photospheric absorption in our models with 
relatively weak winds. 

To estimate the precision and validity of our relations, 
we have rederived the mass loss rates for the 13 stars from 



2 Provided the wind velocity structures are identical (~ same 
/?) this effect should become more important for weaker winds. 
This trend is indeed apparent in Fig. 18b. 



the sample of LL93 having a net observed Ha emission 3 . 
The mass loss rate was directly determined from the ob- 
served Ha equivalent width using Eq. 1 without any pho- 
tospheric or other correction. The observational mass loss 
rates are taken from the new rederivation by Puis et 
al. (1995, their Table 11, column 3), based on the equiv- 
alent widths compiled by LL93. For the individual stars 
the result is plotted in Fig. 19b, the mean relation being 
indicated by the solid line. The same is also shown us- 
ing the stellar parameters from LL93 and Eq. 3 (squares; 
mean relation: dashed). From this figure we see that our 
mean relations yield Ha mass loss rates which are on the 
average 0.1 to 0.3 dex lower than the values derived by 
Puis et al. This is a direct consequence of the differences 
illustrated above. 

5.1.2. The infrared spectrum: Paa & Bra 

The IR recombination lines from the Paschen and Brack- 
ett series show a similar behaviour as Ha (cf. Fig. 16). 
With the advent of new IR observations it may be of in- 
terest to determine the mass loss rate of OB stars from 
the strongest emission lines, namely Paa 1.88 //m and Bra 
4.05 lira. In Fig. 20 we have used predicted strengths to 
plot the mass loss rate as a function of equivalent width, 

3 The smallest equivalent width is of 0.42 A emission. 
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Fig. 20. Mass loss rate as a function of the logarithm of the predicted equivalent width (in A) for Paa (left) and Bra (right). 
Same symbols as in Fig. 1 with 3 additional models (cf. text). The solid lines show the relations given by Eqs. 4 and 5 respectively 



for the same models as for Ha above (again restricted to 
net emission in the lines). Note that we have not included 
the He I (4-3) blend in Paa. Depending primarily on T e ff 
and the He abundance, this would lead to a small increase 
of the total Paa equivalent width and thereby to a flatter 
relation. Similarly the He II (10-8) transition which blends 
Bra was not considered. 

Figure 20 clearly exhibits the strong dependence of 
Paa and Bra on M. The following mean relations are 
obtained from our model results: 

logM = (0. 511 ±0.102) log VF(Paa)- 5. 907 ±0.152 (4) 

logM = (0. 513 ± 0.095) log VF(Bra) - 6. 169 ±0.188 (5) 

The residual rms is 0.18 and 0.17 dex respectively, which 
is of the same order as for Ha. 

In analogy to the case of Ha discussed above, the ac- 
curacy of the M determination can be improved if, in 
addition, R^, T e ff and are known. For the tempera- 
ture correction ct (Eq. 2) the constants cl and c2 have 
the following values: cl = 2.311 (1.303) and c2 = 1.315 
(0.842) for Paa and Bra respectively. The resulting fit 
values for Eq. 3 are a = 0.703 (0.638) and b = -22.111 
(-22.785). We note that, given the large Paa and Bra emis- 
sion, these lines should be less affected by the uncertainty 
found above for Ha from the comparison with Puis et 
al. the models. 

As a first comparison between determinations based 
from the optical and the IR, we calculate M for £ Pup- 
pis (HD 66811). The observed Ha equivalent width of - 
3.09 A is taken from the compilation of LL93. From the 



line profiles of Kaufl (1993) we measured Bra equivalent 
widths of 74 to 101 A. The derived mass loss rate (in 
Mqyi- 1 ) is logM= -5.37 from Ha, and log M= -5.21 to 
-5.14 from Bra, which shows a reasonable agreement be- 
tween both determinations 4 . From their detailed analysis 
including a line profile fit of Ha, Puis et al. (1995) de- 
rive log M= -5.23, which confirms our value obtained from 
pure equivalent widths measurements. If, in addition, we 
take the stellar parameters of LL93 into account and use 
Eq. 3 with the corresponding constants, we obtain logM= 
-5.39 from Ha, and logM= -5.31 to -5.27 from Bra again 
in good agreement with the observed value. More detailed 
comparisons are beyond the scope of the present study. 

6. Summary and conclusions 

We have presented the spectral evolution of main sequence 
stars with initial masses M; = 40, 60 & 85 Mq, us- 
ing "combined stellar structure and atmosphere models" 
(CoStar). The first Co Star studies of the advanced Wolf- 
Rayet phases are presented by Schaerer (1995ab). CoStar 
models consistently combine stellar structure calculations 
with non-LTE line-blanketed model atmospheres. The 
technical details of CoStar and the evolution tracks used 
in this work are discussed in Paper I of this series. The 
main goal of this second paper is to confront our models 



4 Despite the known He enrichment of £ Puppis the use of 
our fit formulae is justified for this case, since the total Ha 
equivalent width including the He II 6520 blend is expected to 
vary by less than ~ 8%. 



a more 'direct' comparison of stellar evolution calculations 
with observations. 

The presented predictions for the spectral evolution 
cover: (a) line blanketed flux distributions from the EUV 
to the far IR, (b) the integrated ionizing fluxes for the 
hydrogen and helium continua, (c) photometry for eleven 
broad-band filters in the optical and infrared and for four 
filters in the UV, (d) synthetic UV spectra (AA/A up to 
10 3 ) for the IUE wavelength range from 1000 to 2200 A, 
including the dominant wind and metal line features, and 
finally (e) the optical and infrared hydrogen and helium 
line spectrum. 

We have performed these calculations at several points 
along the main sequence for three different initial masses. 
They therefore present a homogeneous set of results for a 
large domain of stellar parameters of massive stars, which 
will be of value for studies of nebulae encompassing early- 
type stars and for studies of young stellar populations in 
general (e.g. in H II regions or starbursts). 

At this point we will not summarize all our findings, 
but point out three interesting results. The first concerns 
the ionizing spectrum, the second the iron spectrum, and 
the third is a method to derive mass loss rates from optical 
and infrared lines. 

Comparisons of our CoStar ionizing spectra with other 
calculations based on LTE or non-LTE plane parallel 
models clearly show the importance of wind effects on the 
determination of the ionizing flux of O and, at least, early 
B stars (cf. Gabler et al. 1992, Najarro et al. 1995). This 
may be of considerable importance for the interpretation 
of spectra of H II regions and related objects. 

Since iron is the main contributor to line blanketing 
and to the radiative force in the investigated stars, we 
have paid special interest to diagnostics related to these 
two fundamental problems. The broad band UV colors are 
found to be quite sensitive to line blanketing effects. A 
comparison of these colors with spectrophotometric IUE 
observations shows an encouraging agreement. A number 
of metal line features in the spectral range mostly covered 
by the IUE short wavelength camera are used to investi- 
gate the ionization structure of iron. The strongest fea- 
ture in the spectrum is due to a large number of Fe III 
lines at A ~ 1920A. The predictions of this feature, which 
shows a strong temperature dependence, are in excellent 
agreement with observations. They may even be used as 
a temperature indicator for late O and early B giants and 
supergiants (18 T e ff 33 kK). Other iron features, such 
as the relatively weak blend of Fe V lines around 1453 
A, do not provide such strong diagnostic tools. The Fe V 
1453 feature even points to a somewhat underestimated 
iron ionization in our models. In general however, we con- 
clude that the predictions of various investigated spectral 
features are in reasonable to satisfactory agreement with 
observations. 



to derive simple fit formulae, which allow the determina- 
tion of the mass loss rate from the equivalent width for 
emission of each of these three lines. A test of this method 
using the 04I(n)f star ( Puppis indicates that this method 
yields reliable results. 
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